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. . Abstract 
ON 

The solar contribution to global mean air surface temperature change is analyzed by using an empirical bi-scale climate model 
characterized by both fast and slow characteristic time responses to solar forcing: T] = 0.4 ±0.1 yr, and T2 = 8 + 2 yr or 
T2 = 12 + 3 yr. Since 1980 the solar contribution to climate change is uncertain because of the severe uncertainty of the total solar 
O jrradiance satellite composites. The sun may have caused from a slight cooling, if PMOD TSI composite is used, to a significant 
warming (up to 65% of the total observed warming) if ACRIM, or other TSI composites are used. The model is calibrated only 
on the empirical 11 -year solar cycle signature on the instrumental global surface temperature since 1980. The model reconstructs 
the major temperature patterns covering 400 years of solar induced temperature changes, as shown in recent paleoclimate global 
CN [temperature records. 
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O 1. Introduction 

^ : 

Estimating the solar contribution to global mean air surface 
temperature change is fundamental for evaluating the anthro- 
. ^ pogenic contribution to climate change. This is regarded as one 
^^pf the most important issues of our time. While some theoret- 
(— I ical climate model studies [Hegerl et al., 2007; Hansen et ai, 
Q_|2007; IPCC, 2007] indicate that the solar variability has little 
effect on cUmate (these studies estimate that less than 10% of 
I the global warming observed since 1900 is due to the Sun), sev- 
eral empirical studies suggest that large climatic variations are 
On well synchronized with solar variations and, therefore, climate 
is quite sensitive to solar changes [Eddy, 1976; Hoyt and Schat- 
ten, 1997; White et ai, 1997; van Loon and Labitzke, 2000; 
Douglass and Clader, 2002; Kirkby, 2007; Scafetta and West, 
2005, 2006, 2007, 2008; Shaviv, 2008; Eichler et al, 2009; 
Soon, 2009]. 

I Theoretical studies rely on climate models. Two alternative 
approaches are commonly used: energy balance models (EBM) 
.5^ [for example: Crowley, 2000; Foukal et al, 2004] and general 
circulation models (GCM) [for example, Hansen et ai, 2007]. 
5h These models are based on the idea that climate is forced by 
solar variations, volcano activity, aerosols and several green- 
house gases (CO2, CH4, etc). These forcings are theoretically 
evaluated and used as inputs of the models. The climate sen- 
sitivities to the forcing is estimated according to the known 
physics. This known physics is implemented in the models. 
The models contain a certain number of climate mechanisms 
such as water vapor feedback, cloud formation, energy transfer, 
etc. The major problem with this approach is that the physics 
implemented within the models may be severely incomplete. 
Specifically, some key variables such as the climate sensitivity 
to CO2 changes is severely uncertain. 

For example, according the IPCC [2007] a doubling of CO2 



O 



may induce a temperature increase from 1.5 K to 4.5 K, and 
more. This large uncertainty is mostly due to the current poor 
understanding and modeling of water vapor and cloud forma- 
tion feedbacks which can have large effects on climate [Kirkby, 
2007; Shaviv, 2008]. Indeed, significant discrepancies between 
climate model predictions and data are observed [Douglass et 
ai, 2007; Lean and Rind, 2008], and several climate mecha- 
nisms are still poorly understood, as reported by numerous sci- 
entific papers [Idso and Singer, 2009]. 

An alternative approach is based on empirical multilinear re- 
gression models. It is assumed that not all physics is known or 
implemented in the models. The forcings are used as inputs of 
EBMs whose outputs are not the actual temperature signatures 
generated by the various forcings but waveform functions that 
are assumed to be proportional to such signatures. The temper- 
ature is supposed to be a linear superposition of these rescaled 
output waveforms and linear amplification coefficients are eval- 
uated by means of a multilinear regression analysis of a given 
temperature record. Thus, it is assumed that: 



AT(t) = yaFSF{t) + N(t), 



F 



(1) 



where: the regression coefficients, af, are the linear amplifi- 
cation coefficient associated to a given forcing F; S fit) is the 
output waveform generated by the chosen EBM forced with a 
given forcing F(t); and N(t) is the residual signal that is inter- 
preted as natural climate variability. The above methodology 
has two major variants according to the particular EBM used to 
generate the waveforms. 

Some authors [North et ai, 2004; Hegerl et ai, 2006, 2007] 
use typical EBMs. The adoption of EBMs is particularly useful 
if the interest focuses on local temperature records, but becomes 
less useful if the interest is in the global average temperature. In 
fact, when the EBM outputs need to be averaged on the entire 



Preprint submitted to J. of Atmospheric and Solar-Terrestrial Physics 



December 22, 2009 



globe an EBM does not perform too much differently from a 
simple low pass RC-like filter with appropriate relaxation time 
responses. The relaxation time response of a thermodynamic 
system is related to the heat capacity of the system itself. For 
example, I found that the EBM used by Crowley [2000], where 
the output is averaged on the entire globe, is approximately sim- 
ulated with a low-pass RC-like filter with characteristic time 
T = 10 yr, as deduced from the data published with Crowley's 
paper. In fact, some other authors [for example, Lockwood, 
2008] use low-pass RC-like filters with a specific characteristic 
time response for each forcing. 

On the contrary, other authors [Douglass and Clader, 2002; 
Gleisner and Thejll, 2003; Lean and Rind, 2008] do not use 
traditional EBMs. These authors just assume that the output 
waveform functions coincide with the corresponding forcing 
functions with some time-lag shift. Thus, these authors use Eq. 
[I]with5F(f) = F{t-TF). 

The results of these multilinear regression model studies are 
quite interesting, also because they differ significantly from 
each other. Hegerl et al. [2007] found a large variability of the 
climate sensitivity to the total solar irradiance (TSI) changes 
depending on the paleoclimate temperature records that they 
used. In some case these authors even found negative values 
of the climate sensitivity to TSI changes which is evidently not 
physical because it would imply that global climate cools when 
TSI increases and warms when TSI decreases. Probably, the 
significant uncertainty present in the paleoclimate temperature 
reconstructions and in the forcing functions is responsible for 
these ambiguous results. These results show that the multilin- 
ear regression analysis methodology is inefficient when applied 
to long and uncertain records. 

Lockwood [2008] applied a nonlinear multivariate fit with 
several parameters on a three decades surface temperature 
record and found that the surface climate signature associated to 
the 1 1 -year solar cycle has a peak-to-trough amplitude of about 
0.05 K. On the contrary, Tung and Camp [2008] using similar 
data found a peak-to-trough solar signature amplitude of about 
0.2 K. Douglass and Clader [2002], Gleisner and Thejll [2003], 
Lean and Rind [2008] and several other studies [White et al, 
1997; Scafetta and West, 2005] found that the surface climate 
signature associated to the 11 -year solar cycle has a peak-to- 
trough amplitude of about 0. 1 K. Indeed, this O.IK solar cycle 
signature in the global surface temperature appears to be the 
most common result among the empirical studies [IPCC 2007, 
see page 674 for details], in particular since 1980. Herein, I 
will refer to it as the empirical estimate of the 11 -year solar 
cycle signature on global surface temperature since 1980. 

Indeed, it is relatively easy to find this signature. Figure 1 
shows the original global surface temperature [Brohan et al, 
2006] (curve 'b'), and the volcano (curve 'c') and the ENSO 
(curve 'd') temperature signatures, as recently estimated by 
Lockwood's model [2008]. The curve 'a' in the figure shows 
the temperature detrended of the volcano and of the detrended 
ENSO signature components. The detrended ENSO signature 
component is obtained by detrending the ENSO signature of 
its four year moving average smooth curve, which is shown in 
the figure in the solid thick curve 'd' . This operation does not 
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Figure 1 : Temperature components. The curve (b) is tiie original global surface 
temperature [Brohan et al, 2006]. The curve (c) is the volcano signature on 
the temperature as estimated by Lockwood [2008]. The thin curve (d) is the 
ENSO signature on the temperature as estimated by Lockwood [2008]; the thick 
curve is a four year moving average of the thin curve. The thin curve (a) is 
the surface temperature minus the volcano and ENSO signatures plus the thick 
smooth curve in (d); the thick smooth curve in (a) is a four year moving average 
of the thin curve (a). The curves are dislocated at 0.5 K intervals for visual 
convenience. The "*" symbols indicate the position of the TSI maxima. 

change the final results drastically but it is done because the 
ENSO signature may be capturing part of the solar decadal sig- 
nature on climate, so this smooth component is put back in the 
data before a comparison with the solar record is studied. Also 
Lockwook's residual signal may still contain a solar signature: 
therefore, it is kept in the data to avoid an inappropriate filter- 
ing. 

The filtered temperature signal (curve 'a' in Figure 1) shows 
a clear decadal oscillation with a peak-to-trough amplitude of 
at least 0.1 K, which is in phase with the solar cycles. The "*" 
symbols in the figure indicate the position of the 1 1-year solar 
cycle maxima and, on average, there is a lag-time of about 1 
year between the solar maxima and the maxima of the smooth 
curve 'a', which fits the prediction of some EBMs [see figure 
lb in Northef a/., 2004]. 

The peak to trough empirical amplitude regarding the 11- 
year solar cycle signature on global surface temperature is not 
reproduced by traditional GCM and EBM estimates. North et 
al [2004] used five different EBMs and found that the climate 
signature associated to the 11 -year solar cycle is, on average, 
twice than the theoretical predictions (see their figures 1 and 4). 
The climate models used by Crowley et al [2000], Foukal et 
al [2004] and Hansen et al [2007] predict an even lower solar 
signature on climate with a peak to trough amplitude of about 
0.02-0.04 K. It is reasonable to think that current climate mod- 
els are missing important climate mechanisms that amplify the 
solar signature on climate, also by a large factor [Shaviv, 2008]. 
In fact, these models assume that the sun can alter climate only 
by means of direct TSI forcing while there are strong evidences 
that variation of direct UV radiation and cosmic rays, which af- 
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feet eloud formation and ehange the albedo, can play a major 
role in climate change [Pap et al, 2004; Kirkby, 2007]. Thus, 
there are both empirical and theoretical reasons to believe that 
traditional climate models cannot faithfully reconstruct the so- 
lar signature on cUmate and are significantly underestimating 
it. 

The alternative approach that is based on multilinear regres- 
sion reconstruction of cUmate has also some serious short- 
comings. Multilinear regression analysis is very sensitive to 
the shape of the temperature function and to the shape of the 
functions used as constructors. Thus, uncertainties in the data 
and/or in the models used to construct the waveform com- 
ponents yield suspicious regression coefficients, as Hegerl et 
al. [2007] found. Moreover, multilinear regression analysis 
is based on the assumption of linearity and independence of 
the waveforms. These conditions are not sufficiently satisfied 
particularly for long sequences because, for example, the solar 
forcing, by altering climate, indirectly alters GHG concentra- 
tions too. Thus, GHG forcing cannot be considered indepen- 
dent from solar forcing, as all studies adopting this methodol- 
ogy as well as most EBMs and GCMs have assumed. Con- 
sequently these models underestimate the decadal and secular 
solar contribution to climate change by physical construction. 

Finally, while analyzing shorter time series may reduce some 
uncertainties, when long sequences are analyzed it is funda- 
mental to have a physically accurate model with the correct 
characteristic time responses. We are left with some arbitrari- 
ness: Should we use a model with a characteristic time response 
with a decadal scale, as the EBMs used by Crowley et al. [2000] 
and North et al. [2004] assume, or a model that uses as wave- 
forms the forcing functions with short time-lag shift, as Dou- 
glass and Clader [2002], Gleisner and Thejll [2003] and Lean 
and Rind [2008] do? There is a significant difference between 
the two approaches. In fact, EBMs predict that the time-lag 
shift and the climate sensitivity to a cyclical forcing greatly in- 
crease with decreasing frequency [Wigley, 1988] and, therefore, 
cannot be kept constant at all temporal scales. 

The models used by Douglass and Clader [2002], Gleisner 
and Thejll [2003], Lean and Rind [2008] would be severely 
misleading when applied to multidecadal and secular sequences 
because the amplitude of the low frequency component of the 
solar signature on climate would be severely underestimated. 
The multilinear regression approach with fixed time-lag shifts 
is more appropriate if applied to relatively short time sequences 
(a few decades) and if the study is limited to evaluate the 11- 
year solar cycle signature on cUmate. In fact, if the forcing 
function is characterized by a unique frequency using the latter 
method would be approximately equivalent to using a low-pass 
RC-like filter model or an EBM because a sinus-like function is 
transformed into a fixed time-lag shifted sinus-like function by 
a low-pass RC-like filter or any cUmate model. 

Thus, we have to conclude that both traditional cUmate mod- 
els and multilinear regression analysis models are not com- 
pletely satisfactory. The peak to trough empirical amplitude 
of 0.1''C regarding the 11-year solar cycle signature on global 
surface temperature seems to be sufficiently robust because it 
was obtained with multiple analysis methods by several au- 



thors. Moreover, while climate responds to a given forcing with 
given time responses, the cUmate characteristic time responses 
too should be empirically measured, and not just theoretically 
deduced as in all above studies. 

The climate characteristic time responses were empirically 
measured by Scafetta [2008al, and confirmed by Schwartz 
[2008], by studying the autocorrelation properties of the global 
surface temperature record, assuming that the memory of its 
fluctuations is described by autoregressive models. The result 
is consistent with the fluctuation dissipation theorem that relies 
on the assumption that the response of a system in thermody- 
namic equilibrium to a small applied force is the same as its 
response to a spontaneous fluctuation. 

It was found that cUmate is characterized by two major char- 
acteristic time responses: one short with a time scale of a few 
months (ti - 0.4 + 0.1 yr); and one long with a time scale that 
may be as short as ti = 8 ± 2 yr or, by taking into account sta- 
tistical biases due to the shortness of the available temperature 
record, as long as T2 = 12 ± 3 yr. Thus, the climate system 
appears to be made of two superimposed systems characterized 
with a fast and a slow response to the forcings, respectively. For 
example, the atmosphere has a low heat capacity compared to 
the ocean and may be characterized by this short characteris- 
tic time response. Also cloud dynamics is relatively fast. On 
the contrary, the ocean has a large heat capacity that can be 
responsible for a decadal characteristic time response of the cU- 
mate system. Moreover, other climate phenomena such as the 
change of albedo due to the melting of the glacials and natural 
forestation/desertification processes have a decadal time scale. 
Indeed, several studies confirm that the climate system may be 
characterized by both long and short response times [see ref- 
erences in Lockwood 2008]. The average surface temperature 
record should be the result of the superposition of at least two 
signals generated by the two kinds of cUmate processes. 

Herein, I propose a model that includes the above two char- 
acteristic time responses and I calibrate it on the empirical es- 
timate of the 11 -year solar cycle signature on global surface 
temperature since 1980, which is herein assumed to be 0.1 K. 
Then, this model is used to predict the solar signature on climate 
since 1600 and this is compared with paleoclimate temperature 
reconstructions. 

2. Total solar irradiance records 

Determining how solar activity has changed on decadal and 
secular scales is necessary to estimate the solar contribution to 
cUmate change. Unfortunately, how solar activity has changed 
in time is not known with certainty. 

Direct TSI observations started in 1978 with satellite mea- 
surements. For the period before 1978 only TSI proxy recon- 
structions have been proposed (for example: Hoyt and Schatten 
[1997]; Lean [2000]; Wang et al, [2005] and Krivova et al. 
[2007]). These TSI proxy models significantly differ from each 
other, in particular about the ampUtude of the secular trends. 

Unfortunately, TSI sateUite composites since 1978 are not 
certain either. Two major composites have been proposed: the 
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Figure 2: Three alternative TSI satellite composites. The composite [B] and 
[C] are shifted by -5W/m^ and -lOW/m^, respectively, from the 'native scale' 
for visual convenience, (units of watts/meter^ at 1 A.U.) 

PMOD TSI composite [Frohlich and Lean, 1998; Frohlich, 
2004; Frohlich, 2006] which shows an almost constant trend 
from 1980 to 2000; and the ACRIM TSI composite [Willson 
and Mordvinov, 2003], which shows an increasing trend during 
the same period. 

GCMs and EBMs adopted by the IPCC [2007] assumed that 
TSI did not change significantly since 1950 and that, conse- 
quently, the sun could not be responsible for the significant 
warming observed since 1975. These estimates are based on 
TSI proxy models such as those prepared by Lean [2000] and 
Wang et al. [2005] which are apparently supported by PMOD 
[Frohlich, 2006]. However, the above TSI proxy models would 
be erroneous if the ACRIM TSI composite more faithfully re- 
produces the TSI behavior during the last decades. 

The ACRIM-PMOD controversy is quite complex and, 
herein, a detailed discussion on this topic is not possible. A 
recent work by Scafetta and Willson [2009] reopened the issue 
by providing a careful analysis of the most recent TSI proxy 
model [Krivova et al., 2007] based on magnetic surface fluxes. 
This has been done by estabhshing that a significant degrada- 
tion of ERBE TSI satellite likely occurred during flie ACRIM- 
gap (1989-1992.5), as the ACRIM team has always claimed. 
Moreover, Scafetta and WiUson invaUdated the specific correc- 
tions to Nimbus7 that the PMOD TSI composite requires and 
confirm the opinion of the original Nimbus7 experimental team 
that no sudden increase of the Nimbus7 sensitivity occurred 
on September 29, 1989 [see Hoyt's statement in Scafetta and 
Willson, 2009]. Finally, Scafetta and WiUson [2009] showed 
that the agreement between PMOD and the proxy reconstruc- 
tion about the absence of a trend between the TSI minima in 
1986 and 1996 is coincidental because a careful comparison be- 
tween the proxy model and the unquestioned satelhte data be- 
fore and after the ACRIM-gap proves that the TSI proxy model 
by Krivova et al. [2007] is missing an upward trend. 

Scafetta and Willson [2009] also showed that during the 
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Figure 3: Three alternative secular TSI records constructed by merging the TSI 
proxy reconstruction by Solanki's team [Krivova et at, 2007] with the three 
alternative TSI satellite composites shown in Figure 1. The composite [B] and 
[C] are shifted by -5W/np- and -lOW/m^, respectively, from the 'native scale' 
for visual convenience, (units of watts/meter^ at 1 A.U.) 



ACRIM-gap Nimbus7 may have increased its sensitivity by 
about Q.3W/m^. While this error is well below that hypothe- 
sized by Frohlich (0.86W/m^) to construct PMOD TSI compos- 
ite, this discrepancy is compatible with the known uncertainty 
of Nimbus7 whose sensors, as well as those of ERBE, were not 
able to make accurate self-caUbrations. A direct comparison of 
the local trends between Nimbus7 and ACRIM 1 does show dis- 
crepancies that can be as large as +0.3 W/m^ [Scafetta, 2009]. 

Consequently, herein I assume three TSI composites 
[Scafetta, 2009] that approximately cover the entire range of 
possible TSI satellite composites. The composites are con- 
structed by taking into account overlapping regions and con- 
tinuity at the merging dates. The composite [A] assumes that 
during the ACRIM-gap the Nimbus7 record is accurate; the 
composite [C] assumes that during the ACRIM-gap the ERBE 
record is accurate; the composite [B] is just the arithmetic av- 
erage between [A] and [C] . Note that the ACRIM composite is 
approximately between the composite [A] and [B], while the 
PMOD composite is approximately reproduced by the com- 
posite [C]. All configurations between the composites [A] and 
[C] are ideally possible according to the published TSI satel- 
hte records. As Figure 2 shows, the possible range of TSI 
difference between the two solar minima in 1986 and 1996 is 
Al ^ 0.3 + 0.4W/m2. 

A secular TSI record can be obtained, even if imperfectly, by 
merging a TSI secular proxy reconstruction (here I use the most 
recent TSI proxy reconstruction by Solanki's team [Krivova et 
al. , 2007]) with the three TSI sateUite composites shown in Fig- 
ure 2. The 1980-1990 TSI mean values are used for this merg- 
ing. The three composed TSI records are shown in Figure 3 and 
indicated wifli [A], [B] and [C]. 
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3. An empirical climate model with short and long charac- 
teristic time responses 

As explained in the Introduction, climate appears to be char- 
acterized by at least two major characteristic time constants 
by both theoretical and empirical findings [Scafetta, 2008a; 
Schwartz, 2008]: a short characteristic time constant, ti = 
0.4 + 0.1 yr, and a decadal one from T2 = 8 + 2 yr to T2 = 12 + 3 
yr Consequently, the signature of a given forcing F{t) on the 
global surface temperature, ATf{t), is the superposition of at 
least two major signals: one generated by the processes with a 
short time response and the other generated by those processes 
with a long time response to a forcing. In the case of the solar 
forcing we should write: 



ATsit) = ATis{t) + AT2sit), 



where 

dATis(t) 
dt 

dAT2s(t) 
dt 



kisAI(t)-ATis(t) 
k2sAI(t)-AT2s(t) 

TlS 



(2) 

(3) 
(4) 



The parameters ts\ and tsi are the short and long climate char- 
acteristic time responses, and k\s and k2s are the average equi- 
librium climate sensitivities of the two kind of processes. With 
appropriate parameters, the above model simulates the perfor- 
mance of any energy balance models, as those used by North et 
al. [1981] and Crowley ef a/. [2000]. 

The above equations assume that the TSI record is used as a 
proxy for the overall climate sensitivity to solar changes. Thus, 
the parameters kis and k2s do not have the meaning of climate 
sensitivity to TSI variation as assumed in the EBMs and GCMs, 
but they have a meaning of climate sensitivity to "total solar 
activity" variation in TSI units. The difference is fundamental 
because climate is also altered by solar changes alternative to 
TSI changes. The empirical methodology I propose to evaluate 
kis ™d k2s is summarized in four steps: 

1) We are interested in evaluating the solar signal on climate. 
This is determined by Eq. |2]which depends on four parameters, 
kis, TlS, k2s and T2s- 

2) For the fast process I set tis - 0.4 yr as empirically mea- 
sured [Scafetta, 2008a], and kis = 0.053 KlWm-^. The justifi- 
cation for the latter value is that the climate system has no suf- 
ficient time to physically, chemically and biologically evolve, 
that is, I assume that on this short time scale the chemical com- 
position of the air, the ice cover, the forest cover and all other 
major physical climate variables do not change much. In this 
situation the climate sensitivity can be approximately calculated 



dT 
'dl 



T 
41 



0.053 K/Wm" 



(5) 



where the average Earth surface temperature is T - 289 K and 
the average TSI is / = 1365 W/irP'. The above value is easily 
calculated by differentiating the energy balance law. 



where cr is the Stefan-Boltzmann constant. The albedo "fl^" as 
well as any additional climate function "gr which is needed to 
correct the black-body Stefan-Boltzmann law to make it com- 
patible with the surface Earth's climate system, are assumed to 
be constant on very short time scales. The function "gr" would 
include all feedback mechanisms, emissivity, dispersion of en- 
ergy in the ocean, etc. Thus, for short time scales it is assumed 
that the function h{T) - gril - aj) - const. Note that also 
Lockwood [2008, table 1] found that a climate sensitivity of 
k = 0.052 K/Wm.-^ to solar chan ges is associated to a monthly 
characteristic time response. 

3) The long characteristic time constant of climate is set to 
be as short as T52 = 8 yr or, alternatively, as long as - 12 
yr, as Scafetta [2008a] and Schwartz [2008] have found. 

4) The last free parameter k2s is determined by imposing that 
the peak-to-trough amplitude of the global climate response to 
the 11-year solar cycle is about 0.1 K near the surface after 
1980, as reported by the IPCC [2007, page 674 for details], as 
found by several empirical studies discussed in the Introduction 
and as shown in Figure 1 . 

A reader should be careful here because s/he might mistak- 
enly believe that I am not taking into account the other climate 
forcings such as volcano and GHG forcings. The four parame- 
ters {kis , Tis , k2s and ts2) of the model have been set by using 
findings that do take into account the other forcings. For ex- 
ample, kis is directly calculated from the energy balance law; 
Tsi and Ts2 are measured from the autocorrelation properties 
of temperature residual after detrending a theoretical effect of 
aU forcings [Scafetta, 2008]; the peak-to-trough amplitude of 
the global climate signature to the 1 1-year solar cycle is found 
to be about O.IK near the surface after removing the volcano, 
GHG+aerosol and ENSO signal [Douglass and Clader, 2002; 
Lean and Rind, 2008]. These values have been estimated with 
the best data available since 1980. Thus, by using such values 
as constraints of the model I am taking into account the other 
forcings too, but implicitly, through the constraints. 

Eq. |6]cannot be used directly because it is not known how the 
climate function h(T) changes as function of the temperature. 
Thus, I use the O.IK solar cycle signature near the surface as a 
natural oscillatory output that enables us to determine the value 
of the climate modulation transfer function to the 1 1 -year solar 
variation, and from this I can determine k2s , as explained below. 

The climate sensitivity to the 11 -year solar cycle is about 
Z(ll) = 0.11 + 0.02 K/Wm-^ [Douglass and Clader, 2002; 
Scafetta and West, 2005], that is, an 1 1-year solar cycle of about 
1 W/rn^ causes a peak to trough cycle of about 0.1 "C on the 
global surface. Thus, Eq. (|2| is solved by the system: 



Z(ll) 
kis 

k2S 



Zi(ll) + Z2(ll) 




(7) 
(8) 

(9) 



4cr 



h(T)I, 



(6) 



where '11' refers to the period of the 11 -year solar cycle, and 
Zi ( 1 1 ) and Z2 ( 1 1 ) are the climate sensitivities to the 1 1 -year so- 



5 





ESS = ESS1 + ESS2 
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Figure 4: An example of ESS curves. ESSl curve shows that the empirical 
solar signature on climate obtained with Eq. (8) with T[ = 0.4 yr and k[s = 
0.053K/Wm-^. ESS2 curve shows that the empirical solar signature on climate 
obtained with Eq. |9) with T2 = 12 yr and k2s = QAIK/Wm-^. ESS (black 
line) shows the superposition of ESSl and ESS2. The model is forces with 
the TSI reconstruction [B] shown in Figure 3. The figure shows temperature 
anomalies relative to the 1895-1905 average. 



lar cycle relative to climate processes with fast and slow char- 
acteristic time responses, respectively. 

The above system can be easily solved by evaluating first 
Zi(ll), then Z2(ll) and, finally, k2s- By using Eq. ^ 
with Ti = 0.4 yr and kis = 0.053 K/Wm-^, it is found: 
Zi(ll) =s 0.051 K/Wm-^ andZ2(ll) ~ 0.059 K/Wm.-^. The 
former value is about 46% of the empirical climate sensitiv- 
ity to the 11-year solar cycle Z(ll) = 0.11 ± 0.02 K/Wm-^. 
Thus, the peak-to-trough amplitude of the global climate re- 
sponse to the 11 -year solar cycle associated to the fast pro- 
cesses is about 0.046 K near the surface against the empiri- 
cal value of 0.1 K. If r2 ~ 8 yr, then k2s ~ 0.28 K/Wm.-^; 
if T2 ~ 12 yr, then k2s ~ 0.41 K/Wm^^. With the above values 
the overall climate sensitivity to solar changes at equilibrium is 
ks ~ kis + k2s ~ 0.33 K/Wm-^ and ks ~ 0.46 K/Wm-^, re- 
spectively. For convenience, the parameters used in the model 
are summarized in Table 1 . 



4. The empirical solar signature on climate 

Figure 4 shows the empirical solar signature (ESS) on cli- 
mate under the assumptions and scenarios mentioned above. 
ESSl and ESS2 are the empirical solar signatures induced by 
the climate processes with fast, r^i, and slow, ts2, character- 
istic time responses, respectively. The curves ESSl, ESS2 and 
their superposition are shown in the particular case in which 
T2 = 12 yr and the TSI record [B] is adopted. 

The curve ESSl shows a very small multidecadal and secu- 
lar variability. According to ESSl, the sun has induced about 
H-0.05/r of warming from 1900 to 2000. Because the global 
warming observed since 1900 has been about 0.8^r, the sun 
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Table 1: The parameters used in the empirical bi-scale climate model Eqs. |2] 
for two alternative values of T2. 



would have caused about 6% of the observed warming through 
those mechanisms that respond quickly to changes in TSI. This 
result does not differ significantly from the IPCC [2007] aver- 
age estimates obtained with the current climate models. 

However, the small secular variability shown in ESSl is in- 
sufficient to explain the preindustrial secular variability of the 
temperature, as estimated by the paleoclimate proxy tempera- 
ture reconstructions [North et ai, 2006]. These paleoclimate 
temperature reconstructions show that the preindustrial tem- 
perature variability (from 1600 to 1900) ranges from 0.05 to 
Q.5K, while ESS 1 in Figure 4 shows a preindustrial temperature 
variability of about 0.Q4K. The small GHG variations observed 
before 1900 and volcano activity could fill the difference only in 
the extreme case that the preindustrial climate is almost stable, 
as in the hockey stick temperature graph by Mann et al. [1999]. 

However, the hockey stick temperature graph is unlikely be- 
cause the Medieval Warm Period and the Little Ice Age, which 
occurred during a maximum and a minimum of solar activity, 
respectively, are supported by numerous historical facts [Hoyt 
and Schatten, 1997] and data from several regions of the Earth 
[Loehle and McCulloch, 2008]. A global warming between the 
periods 1650-1700 and 1850-1900, at least equal to the aver- 
age among the paleoclimate temperature reconstructions and 
more, {Q3K < AT < Q.5K), can be considered more realis- 
tic. Thus, unless most of the paleoclimate proxy temperature 
reconstructions are severely erroneous and the historical evi- 
dences are severely misleading. Figure 4 further stresses that 
the processes generating ESSl solar signature on climate alone 
cannot explain climate change on a secular scale. 

Figure 4 shows that most of the multidecadal and secular 
variability is captured by ESS2 curve. The climate processes 
with a slow time response to TSI changes are necessary to cor- 
rectly interpret the multidecadal and secular climate changes, 
also since 1980. This large sensitivity is mostly due to the fact 
that the climate function h{T), that appears in Eq. |6] is tem- 
perature dependent. On large time scales its temperature de- 
pendency cannot be neglected because it has a large impact on 
climate. 

For example, if a 1 W/ni^ TSI increase causes at equilib- 
rium an increase of the function h(T) by just 0.5% it is easy 
to calculate that the climate sensitivity to TSI change would be 
ks = ST/61 = OAIK/Wm-^. This value is eight times larger 



than k 



IS 



0.053 K/Wm as it was calculated by assuming 
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Figure 5: ESS curves against the filtered global surface temperature record 
shown in the curve (a) of Figure 2. [A] The ESS curves are obtained with ri = 
0.4 yr, T2 = 12 yr and kis = 0.053K/Wm-^ kis = OAIK/Wm-^. [B] The ESS 
curves (blak) are obtained with n = 0.4 yr, r2 = 8 yr and kis = 0.052K/Wm^^ 
k2s = 0.2SK/Wm~^. The model is forces with the TSI reconstruction [A], [B] 
and [C] shown in Figure 3. The figure shows temperature anomalies relative to 
the 1895-1905 average. 



that the function hiT) is constant. This sensitivity is easily cal- 
culated with the following equation: 



To + dT 



To 



h{To + 6T) lo + 61 



h{To) 



(10) 



where Tq - 289 K and /q = 1365 WjirP-. So, even small temper- 
ature dependency of the energy-balance equation parameters, 
such as the albedo and the emissivity, can amphfy the climate 
sensitivity to solar changes by a large factor, as also Shaviv 
[2008] noticed. 

Figure 4 also shows that the decadal variability associated to 
the 11 -year solar cycle derives almost evenly from both kind 
of climate processes with fast and slow characteristic time re- 
sponses, as reproduced in the ESSl and ESS2 curves. 

Figures 5A and 5B show the ESS curves obtained by super- 
imposing ESSl and ESS2 curves under different assumptions 



and scenarios against the global surface temperature record 
since 1950 detrended of the volcano and ENSO signatures (see 
curve 'a' in Figure 1). Figures 5A and 5B use T2 = 12 yr and 
T2 = 8 yr, respectively. According to the figures there would 
be no significant solar induced warming since 1950 if the TSI 
record fC] is used; in this case the trend would be slightly nega- 
tive since 1980, as also Lockwood [2008] found using PMOD. 
Instead, the sun can induce as much as 65% of the observed 
warming if the TSI record [A] is used. On average the sun may 
have induced a significant warming since 1950 and since 1980, 
respectively. Figures 5A and 5B show also the good correspon- 
dence between the 1 1 -year solar cycle signatures and the cycles 
observed in the four year moving average smooth curve of the 
filtered global surface temperature record. 

Finally, Figure 6 shows a four century comparison between 
two ESS curves and a paleoclimate temperature reconstruction 
[Moberg et al, 2005] from 1600 to 1850 and the global surface 
temperature record since 1850. The two ESS curves shown in 
the figure are those with the highest secular variability (curve 
#1: T2 = 12 yr and with TSI [A]), and with the lowest secu- 
lar variability (curve #2: T2 = 8 yr and with TSI [C]), respec- 
tively. The figure shows that the ESS signals reproduce quite 
well the cooling and warming patterns observed in the temper- 
ature record for four centuries, in particular with the ESS curve 
#1. Since 1980 there is a wide fork whose extremes depend on 
which TSI satelUte composite is used. If the TSI composite [C] 
is used, the sun would have caused just a slight cooling. If the 
TSI composite [A] is used, the sun would have caused a signif- 
icant warming. Thus, on average it is not unlikely that the sun 
has induced a significant warming since 1980 as it was inferred 
by Scafetta and Willson [2009]. 

5. Conclusion 

Herein I have analyzed the solar contribution to global mean 
air surface temperature change. A comprehensive interpretation 
of multiple scientific findings indicates that the contribution of 
solar variabihty to climate change is significant and that the 
temperature trend since 1980 can be large and upward. How- 
ever, to correctly quantify the solar contribution to the recent 
global warming it is necessary to determine the correct TSI be- 
havior since 1980. Unfortunately, this cannot be done with cer- 
tainty yet. The PMOD TSI composite, which has been used 
by the IPCC and most climate modelers, has been found to be 
based on arbitrary and questionable assumptions [Scafetta and 
Willson, 2009]. Thus, it cannot be excluded that TSI increased 
from 1980 to 2000 as claimed by the ACRIM scientific team. 

The IPCC [2007] claim that the solar contribution to climate 
change since 1950 is negligible may be based on wrong so- 
lar data in addition to the fact that the EBMs and GCMs there 
used are missing or poorly modeling several climate mecha- 
nisms that would significantly amplify the solar effect on cU- 
mate. When taken into account the entire range of possible TSI 
satellite composite since 1980, the solar contribution to chmate 
change ranges from a sUght cooling to a significant warming, 
which can be as large as 65% of the total observed global warm- 
ing. 
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Figure 6: ESS curves against a paleoclimate temperature reconstruction 
[Moberg et al, 2005] from 1600 to 1850 (thin grey line) and global sur- 
face temperature record since 1950 [Brohan et al, 2006] (thick black line). 
The ESS (thick grey) curve is obtained with ri = 0.4 yr, tj = 12 yr and 
kis = 0.053A:/lVm-2 kis = OMKIWm-^ with the model forced with the TSI 
reconstruction [A]. The ESS (thin black) curve is obtained with ri = 0.4 yr, 
T2 = 8 yr and kis = 0.053A:/H'm"^ hs = O.ISK/Wm'^ with the model forced 
with the TSI reconstruction [C]. 

The above wide range strongly contrasts with some recent 
estimates such as those found by Loclcwood [2008], who cal- 
culated that the solar contribution to global warming is negligi- 
ble since 1980: the sun could have caused from a -3.6% using 
PMOD to a H-3. 1 % using ACRIM. In fact, Lockwood's model is 
approximately reproduced by the ESSl curve that refers to the 
solar signature on climate as produced only by those processes 
characterized with a short time response to a forcing. Indeed, 
the characteristic time constants that Lockwood found with his 
complicated nonlinear multiregression analysis are all smaller 
than one year (see his table 1) and the climate sensitivity to 
TSI that he found is essentially equal to my kis ' Likely, Lock- 
wood's model was unable to detect the climate sensitivity to 
solar changes induced by those climate mechanisms that have 
a decadal characteristic time response to solar forcing: mecha- 
nisms that must be present in nature for physical reasons. As 
proven above, these mechanisms are fundamental to properly 
model the decadal and secular trends of the temperature be- 
cause they yield high climate sensitivities to solar changes. 

Analogously, my findings contrast with Lean and Rind 
[2008], who estimated that the sun has caused less than 10% 
of the observed warming since 1900. The model used by Lean 
and Rind, like Lockwood's model, is not appropriate to eval- 
uate the multidecadal solar effect on climate. In fact. Lean 
and Rind do not use any EBM to generate the waveforms they 
use in their regression analysis. These authors assume that the 
temperature is just the linear superposition of the forcing func- 
tions with some fixed time-lags. They also ignore ACRIIVI TSI 
satellite composite. While Lean and Rind's method may be 
sufficiently appropriate for determining the 11 -year solar cy- 
cle signature on the temperature records there used, the same 



method is not appropriate on multidecadal scales because cli- 
mate science predicts that time-lag and the climate sensitivity 
to a forcing is frequency dependent. Consequently, as Lock- 
wood's model. Lean and Rind's model too misses the larger 
sensitivity that the climate system is expected to present to so- 
lar changes at the decadal and secular scales. 

I have shown that the processes with a long time response to 
climate forcing are fundamental to correctly understanding the 
decadal and secular solar effect on climate (see ESS2 curve). 
With simple calculations it is possible to determine that if the 
climate parameters (such as the albedo and the emissivity, etc.) 
change slowly with the temperature, the climate sensitivity to 
solar changes is largely amplified as shown in Eq. [10] 

This finding suggests that the climate system is hypersensi- 
tive to the climate function h(T) and even small errors in mod- 
eling h(T) (for example, in modeling how the albedo, the cloud 
cover, water vapor feedback, the emissivity, etc. respond to 
changes of the temperature on a decadal scale) would yield 
the climate models to fail, even by a large factor, to appro- 
priately determine the solar effect on climate on decadal and 
secular scale. For similar reasons, the models also present a 
very large uncertainty in evaluating the climate sensitivity to 
changes in CO2 atmospheric concentration [Knutti and Hegerl, 
2008]. This large sensitivity of the climate equations to phys- 
ical uncertainty makes the adoption of traditional EBlVIs and 
GCMs quite problematic. 

About the result depicted in Figure 6, the ESS curve has been 
evaluated by calibrating the proposed empirical bi-scale model 
only by using the information deduced: 1) by the instrumen- 
tal temperature and the solar records since 1980 about the 11- 
year solar signature on climate; 2) by the findings by Scafetta 
[2008a] and Schwartz [2008] about the long and short charac- 
teristic time responses of the climate as deduced with autore- 
gressive models. The paleoclimate temperature reconstructions 
were not used to calibrate the model, as done in Scafetta and 
West [2007]. Thus, the finding shown in Figure 6 referring to 
the preindustrial era has also a predictive meaning, and implies 
that climate had a significant preindustrial variability which is 
incompatible with a hockey stick temperature graph. 

Acknowledgment: NS thanks the Army Research Office for 
research support (grant W91 lNF-06-1-0323). 
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